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ABSTRACT 

Recent observations show that planet can reside in close binary systems with 
stellar separation of only ~20 AU. However, planet formation in such close bi- 
nary systems is a challenge to current theory. One of the major theoretical prob- 
lems occurs in the intermediate stage— planetesimals accretion into planetary 
embryos— during which the companion's perturbations can stir up the relative 
velocites(Ay) of planetesimals and thus slow down or even cease their growth. 
Recent studies shown conditions could be even worse for accretion if the gas disk 
evolution was included. However, all previous studies assumed a 2-dimentional 
(2D) disk and a coplanar binary orbit. Extending previous studies by including a 
3D gas disk and an inclined binary orbit with small relative inclination of zb=0.1- 
5 deg, we numerically investigate the conditions for planetesimal accretion at 1-2 
AU, an extension of the habitable zone(~l-1.3 AU), around a Centauri A in 
this paper. Inclusion of the binary inclination leads to: (1) differential orbital 
phasing is realized in the 3D space, and thus different- sized bodies are separated 
from each other; (2) total impact rate is lower, and impacts mainly occur between 
similar-sized bodies; (3) accretion is more favored, but the balance between accre- 
tion and erosion remains uncertain, and the "possible accretion region" extends 
up to 2AU when assuming an optimistic (critical specific energy that leads 
to catastrophic fragmentation); and (4) impact velocities (A^) are significantly 
reduced but still much larger than their escape velocities, which infers that plan- 
etesimals grow by means of type II runaway mode. As a conclusion, inclusion of a 
small binary inclination is a promising mechanism that favors accretion, opening 
a possibility that planet formation in close binary systems can go through the 
difficult stage of planetesimals accretion into planetary embryos. 

Subject headings: methods: numerical — planetary systems: formation 
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1. INTRODUCTION 

As of today, more than 300 extrasolar planets have been found around main sequence 
stars, and about 33 of them are known to reside in binary systems(Eggenberger et aL 2004; 
Konachi 2005; Mugrauer et aL 2005). As solar- type stars usually born in binary or multiple 
systems(Duquennoy & Mayor 1991; Mathieu et al. 2000; Duchene et al. 2007), more and 
more planets are expected in binary systems with the improvement of the detecting tech- 
niques. In most detected planet-bearing binary systems, planet resides around one star with 
the other companion acting as a perturber on an outer orbit (so called S-Type orbit). The 
orbital separations of the binary are found in a wide range from thousands to only a few 
tens of AU. Three systems— GJ86 (Queloz et al. 2000), 7 Cephei (Hatzes et al. 2003), and 
HD 41004 (Zucker et al. 2004)— have stellar semimajor axes of only ~20 AU. Moreover, the 
star HD 188753A, which might host a Jupiter-mass planet (supported by Konachi 2005, but 
challengaed by Eggenberger et al. 2007) on a 3.35 day S-Type orbit, has a companion with a 
close periastron distance of only ~6 AU. The discoveries of these close binary systems chal- 
lenge the current planet formation theory with the question that whether planet formation 
can go through without being prevented by the companion's strong perturbations in such 
close binary systems. 

Current standard theory— core accretion model— describes the process of planet forma- 
tion in 3 phases (Lissauer 1993; Papaloizou & Terquem 2006; Armitage 2007): (1) formation 
of km-sized planetesimals(10^^-10^° g) from grain coargulation (Weidenschilling & Cuzzi 
1993) or from gravitational fragmentation of a dense particle sub-disk near the midplane 
of the protoplanetary disk (Goldreich & Ward 1973), on timescales of the order of 10^ yr; 
(2) accretion of planetesimals into planetary embryos (10^^-10^^ g. Mercury- to Mars- size) 
through a phase of runaway and oligarchic growth on a timescale of the order of 10^-10^ yr 
(Greenberg et al. 1978; Wetherill & Stewart 1989; Barge & Pellat 1993; Kokubo & Ida 1996, 
1998, 2000; Rafikov 2003, 2004); and (3) giant impacts between embryos, producing full-sized 
(lO^'^'-lO^^ g) terrestrial planets in about 10''-10^ yr (Chambers & Wetherill 1998; Levison & 
Agnor 2003; Kokubo et al. 2006). For the first stage (grain-planetesimal phase), a recent 
study by Pascucci et al. (2008) suggests that companion stars at projected separation above 
10 AU do not appreciably affect the degree of grain growth. For the later stage(embryo- 
planet phase), simulations show binary with periastron of above 10 AU can form terrestrial 
planets over the entire range of orbits allowed for single stars(Quintana & Lissauer 2007; 
Quintana et al. 2007). Stellar companions even hasten the accretion process by stirring 
up the planetary embryos to higher eccentricities and inclinations which thus reduce their 
orbit-crossing timescales(Quintana 2004; Quintana & Lissauer 2006). For the intermediate 
stage(planetesimal-embryo), however, planetesimal accretion in close binary systems is not 
as easy as that around single stars. 
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The influence of a companion star on the planetesimal accretion was first studied without 
including gas drag by Heppenheimer(1974, 1978) and later by Whitmire et al.(1998) . Their 
results implies that the companion star, especially if it is on a close orbit with high eccen- 
tricity, may prevent planetesimal accretion by exciting high relative velocities(Al^) between 
colliding planetesimals. AV is crucial parameter, determining whether accretion or erosion 
dominates. If AV is greater than the planetesimal escape velocity [Vesc ~ 100i?p/(100km) 
ms~^], runaway accretion is prevented and growth is slowed down. Furthermore, if AV^ is 
excited beyond the threshold velocity Kro(a few times of Vesc, depending on the prescription 
for collision, see the Appendix for details) at which erosion dominates accretion, planetes- 
imal growth ceases. Fortunately, in the early stages, planetary disk was rich in gas, which 
could partly damp companion's perturbations toward an equilibrium. Furthermore, Marzari 
& Scholl (2000) found, in this equilibrium state, all the orbits of planetesimals with the same 
size and similar semi-major axes had the configuration of periastron-aligned (so called "or- 
bital alignment"), hence reducing AV to very low values. Applying this mechanism to the 
a Cen AB system, accretion zone of km-sized bodies could extend to 2 AU to the primary 
star {a Cen A). Nevertheless, their work only took into account the AV^ of the equal-sized 
bodies. Detailed study by Thebault et al. (2006) found the orbital alignment simultaneously 
increased AV of different-sized bodies, and it was so efficient that leaded to significant AV 
increase for any departure from the exact equal-sized condition(i?i = R2, where -Ri and R2 
are the radii of the two colliding bodies). Extending their studies by including the effect of a 
fast gas-disk depletion procedure under UVO photoevaporation from the host star or nearby 
stars(Matsuyama et al. 2003; Alexander et al. 2006), Xie and Zhou (2008) found the effect 
of differential orbital alignment could be reduced and all planetesimals eventually converged 
toward the same forced orbits, regardless of size. This mechanism indeed reduces all AV 
to low values, however it is not efficient enough for small bodies(radii less than 5 km) and 
it take a long time(~ 10^-10^ yr) to reach the "low AV^" state, before which conditions are 
still accretion-hostile. Recently, taking into account the differential phasing effect, Thebault 
et al. (2008, 2009) revisited the planetesimal accretion conditions in the a Cen AB system, 
and found the "accretion-safe" area never extended beyond 0.75 AU around Cen A or 0.5 
AU aournd Cen B. Their results also confirmed that the planetesimals-to-embryos phase was 
more affected by the binary environment than the last stages, and it might be one of the 
biggest difficulties of planet formation in close binary systems. 

All previous studies belong to coplanar cases, which have not taken into account the 
inclination of the companion's orbit relative to the planetary disk. According to Hale (1994), 
coplanarity is just an approximation with a deviation of the relative orbital inclination of 2^ < 
10 degrees for solar-type binary systems with separations less than 30-40 AU, and systems 
with separations greater than 40 AU appear to have their orbital planes randomly distributed. 
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Therefore, a small inclination between the binary orbit and the planetary disk(such as zb=1°- 
5°) is more realistic than exact coplanarity(zB = 0) for the binary systems with separations 
of ~ 20 AU. In this paper, assuming a small inclination for the binary orbit, we reinvestigate 
the conditions for planetesimal accretion around a Cen A. Results are compared with the 
previous study Thebault et al. (2008). As presented next, by considering the inclination of 
the companion star, planetesimal accretion is much more favored than in the coplanar case, 
thus providing a possible solution for the problem of planet formation in close binaries. 

In §2, we describe our numerical model and initial conditions. In §3, we first compare 
two nominal cases: = and = 1 deg, then present the results of other cases. Some 
further discussions are in §4. Finally, we summarize in §5. 



2. NUMERICAL MODEL AND INITIAL CONDITIONS 

2.1. Gas Disk Model 

Gas drag is modelled by assuming a non-evolving axisymmetric gas disk with constant 
circular streamlines. Following Weidenschilling & Davis (1985), the drag force is expressed 

as 

F = -Kvw, (1) 

where F is the force per unit mass, v is the relative velocity between the planetesimal and 
gas, V is the velocity modulus, and K is the drag parameter, defined as 

i^=^^. (2) 

Here Pg is the local gas density, and Pp and are the planetesimal density and radius, 
respectively. Cd is a dimensionless coefficient related to the shape of the body (~ 0.4 for 
spherical bodies). 

As we consider the case of is 7^ 0, both the distribution and motion of gas should be 
modeled in the 3-Dimentional space. According to Takeuchi & Lin (2002), in cylindrical 
coordinates(r, z), we adopt the gas density distribution as 

2 

T Z 

Pg{r,z) = pgo{—Yexp{-—), (3) 

and the gas rotation as 

fi,(r, z) = + li^rip + + (4) 



- 5 - 



where hg{r) = hQ{^Y'^^^^^'^ is the scale height of gas disk, ^lk,mid is the Keplerian rotation 
on the mid-plane, and the subscript "0" means the value at lAU. In this paper, we take the 
Minimum Mass of Solar Nebula (hereafter MMSN, Hayashi 1981) as the nominal gas disk, 
where p = -2.75, q = -0.5, pgo = 1.4 x 10"^gcm"^, and ho = 4.7 x 10"^ AU. 

As stated in many previous studies(Scholl et al. 2007; Thebault et al. 2006), this 
axisymmetric gas disk model is a crude simplification, and in reality, perturbations on the 
gas disk by the companion should be included (Artymowicz & Lubow 1994; Kely & Nelson 
2007). A full description of the dynamical behavior of planetesimals in such a "real" disk 
requires hydro-code modelling of the gas in addition to N-body type models for planetesimals. 
However, there still exist some uncertainties even in a "real" disk model: how to choose a 
numerical wave damping procedure(Paardekooper et al. 2008). On the other hand, such 
an all-encompassing gas-plus-planetesimals modelling costs much in calculation and goes 
beyond the scope of our study in this paper. Fortunately, the simplified axisymmetric gas 
disk model, first of all, is valid on the average; second, it costs very little CPU time and is 
thus suited for detailed parameter-space explorations. Furthermore, as we discuss later, the 
behavior of should probably be similar in both the simplified and full gas models, as 
long as we have considered the binary inclination. 



2.2. Initial Set Up 

Table 1 lists most parameters for a Cen AB system, planetesimals and the gas disk. 
lO'' planetesimals with random eccentricities of - 10"^, random inclinations of - 5 x 10"^ 
rad, random semimajor axes of 1 - 2 AU, and random radii of 1 - 10 km, orbit around 
the primary star, a Cen A. For the nominal gas disk(lMMSN), we explore four cases with 
different binary inclinations from to 5 deg; for the nominal binary inclination (z^ = 1 
deg), we explore three cases with different gas densities(0.1, 1, and 10 MMSN). The inflated 
radius for collision search is set as 5 x 10~^ {vcoi/ AUy^'^ AU, where Vcoi is the distance from a 
collision to the primary. This choice gives a precision of ~ 1.5 m.s~^ in every impact velocity 
estimate. All the cases evolve 10^ yr by using the fourth-order Hermite integrator (Kokubo 
et al. 1998), including the gas drag force and companion's perturbations. 
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3. RESULTS 

3.1. Coplanar Case vs Inclined Case 

In order to clearly see and fully understand the role of the binary inclination, we present 
a detail comparison between the nominal coplanar and inclined cases: ib = and = 1 
deg. Both cases have a nominal gas disk of IMMSN. 

3.1.1. Orbital distribution 

In Figure 1, we plot the evolutions of planetesimals' distribution in the X-Z coordinate 
plane for both = and = 1 deg cases. As can be seen in the = case, at be- 
ginning t = 0, planetesimals with random radii of 1-10 km are randomly mixed together. 
As planetesimals move in the gas disk, they are subject to a gas drag force, which progres- 
sively damps their orbital inclinations. This inclination damping process is more efficient 
for small bodies because gas drag force is inversely proportional to planetesimal radius. At 
the end of the run t = 10^ yr, the average Z value of 1 km bodies is almost one order of 
magnitude lower than those of 10 km ones. Small bodies are crowed in the mid-plane of 
the disks while big ones are diffused on the surface. For the = 1 deg case, although 
with the same initial condition as in the iB = case, the evolution of planetesimals' orbits 
is entirely different. Planetesimals' inclinations are not damped by the gas drag force but 
excited by the companion star (see the second panels on the right side). Under the coupled 
effect of the companion's secular perturbation and gas drag force, different-sized bodies are 
eventually forced to orbits with different orbital inclinations. Because of the size- dependence 
of gas drag force, smaller planetesimals are damped more by gas drag and thus have smaller 
forced inclinations than bigger bodies. As can be seen in the two bottom panels on the right 
side of Figure 1, different- sized bodies are grouped and separated by their radii because of 
the different forced inclinations. In the next two subsections, we will show that this orbital 
grouping significantly changes the distribution of impacts among a swam of planetesimals, 
reduces their impact velocities AV^ and finally favors planetesimal accretion. 

3.1.2. Impact rate distribution 

As found by Thebault et al. (2006), AV is very sensitive to the difference of radii 
sizes(i?i and R2) of the two colliding bodies because of the size-dependence of the orbital 
alignment. In order to understand the evolution and distribution of AV^ among a swarm of 
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planetesimals, we, therefore, should first investigate the impact distribution on the R1-R2 
plane. In Figure 2, we track the evolution of the impact distributions for both ib = and 
is = ^ deg cases. In the first 5 x 10^ yr, the distributions are random and similar in both 
cases because they start with the same random initial condition and orbital alignments are 
not fully completed (see Figure 1). Afterwards, the distributions become entirely different: 
for is = case, impacts occur much more often between bodies of different sizes; for = 1 
deg case, impacts mainly occurs between similar-sized bodies. These two totally different 
impact distributions can be understood as follows: 

i) In the is = case, planetesimal disk is very flattened and can be treat as a 2- 
dimensional (2D) disk. Because of orbital alignment in the 2D disk, a planetesimal can 
collide with another similar-sized one only if their semi-major axes are very close since they 
have the similar forced orbits (means the similar orbital shapes and configurations) and 
radial drifts; in contrast, given the different forced orbits and radial drifts, a planetesimal 
can cross more orbits of other different-sized bodies over a much larger region. 

ii) In the 2^ = 1 deg case, orbital alignment is in 3D space and thus bodies of different 
sizes are separated from each other (see Figure 1), leading to much lower probability of 
impact between two planetesimals of different sizes. 

In addition, a big enhancement of impact rate is observed in the is = case since gas 
drag damps planetesimals' inclinations and thus disk becomes more dense. While in the 
iB = I deg case, impact rate is much lower because planetesimals' inclinations are excited 
and thus the volume density of disk is reduced. 

3.1.3. Accretion vs erosion 

We plot Figure 3 to finally analyze the accretion or erosion conditions for the cases of 
= and = 1 deg. In the two top panels, as expected from the analysis of the last 
subsection, the impacts in = case are much more than those in = 1 deg case. Since 
impacts between different-sized bodies are favored in the = case, the impact velocities 
are increased to relatively high values with average of ~ 200 ms~^; while in the is = ^ deg 
case, impacts mainly occur between similar-sized bodies with low impact velocities of 10 - 
100 ms~^ in average. In the two bottom panels, by comparing each AV^ with two threshold 
velocities: Viow and Vhigh (see the Appendix for detail), we plot the relative weights of the 
possible types of collision outcomes as a function of distance to the primary star (a Cen A). 
For the case of zb = 0, most impacts lead to erosion(red zone). The weight of accretion 
impact (green + blue zone) is only less than 5% or about 10% if the uncertain ones (yellow 
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zone) are included. These results are roughly comparable with the Figure 3 in Thebault et 
al. (2008), which also adopted the assumption of = but a Maxwellian distribution for 
planetesimal size. In contrast, for the case of = 1 deg, accretion is significantly favored: 
at 1.15 AU, center of the habitable zone(~l-1.3 AU, e.g. Barbieri et al. 2002), about 40% 
(80%, if including uncertain ones) collisions lead to accretion. If we define the "possible 
accretion zone" as green+blue+yellow zones >50%, it extends to as far as ~2 AU, beyond 
which companion's perturbations is too strong to reach orbital alignment for planetesimal of 
any size. Note that our "possible accretion zone" here is an optimistic definition, where all 
uncertain ones are taken into account as accretion collision. Actually, the balance between 
accretion and erosion has a large uncertainty (see the yellow region in Fig. 3) because of the 
poor constraint on (critical specific energy that leads to catastrophic fragmentation, see 
the Appendix for detail). 

3.2. Other Cases 

Beside the nominal case with is = 1 deg and gas disk density of IMMSN, we explore 
some other cases with different is and different gas densities. 

3.2.1. Different iB 

Figure 4 shows the results of two cases with the same gas density of 1 MMSN, but 
different binary inclinations: = 0.1 and 5 deg respectively, is = 0-1 deg is marginal case, 
where every planetesimal has a forced inclination in a very small range ~ 0.1 deg and 
thus very close to each other. In such case, bodies of different sizes are not well separated 
from each other and impacts between different-sized bodies are still an important part (see 
the left up panel of Fig. 4). Therefore, the average impact velocity is intermediate [AV ~ 
100 ms~^), leading to, although more accretion-friendly than the case of is = 0, still an 
accretion-hostile condition(red zone exceeds 50% at least). For the larger inclination case of 
is = 5 deg, although the planetesimal vertical excursions might actually exceed the typical 
thickness of the gas disk, the results are similar to the nominal case of = 1 deg, except 
for less impact events because of the larger inclination(see the 3 bottom panels in Fig. 4). 
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3.2.2. Different gas disks 

We also explore two cases with gas densities of 0.1 MMSN and 10 MMSN in Figure 
5. By comparing these two cases with the nominal inclined case in Figure 3, we find: as 
gas density increases from 0.1 to 10 MMSN, the weight of accretion increases in most outer 
region but decreases in the inner region(say the habitable zone ~ 1.0-1.3 AU). The probable 
explanations are: 

i) In the inner region, the companion's perturbations are weak and gas is dense; while 
in the our region, perturbations are stronger and gas density is lower. As a result of these 
conditions, orbital alignment is realized much more easily in the inner than in the outer 
region. 

ii) For 0.1 MMSN gas density, orbital alignment is already realized in the inner region 
but not in outer region. A larger gas density still favors accretion in the outer region, but 
it can no longer improve the inner accretion condition. Furthermore, gas drag is a double- 
edged sword. If orbital alignment is already realized, increasing gas density can result in 
larger different radial drifts between different-sized bodies, which increase AV and finally 
prevents accretion. 

4. DISCUSSION 

4.1. Behaviors In The Full Gas Model 

As mentioned in §2.1, our gas drag model is a simplified one where gas disk is assumed to 
be fully axisymmetric and non-evolving. A full gas model requires hydro-code for modeling 
the gas evolution and N-body type models for planetesimals. Recently, taking the full gas 
model, Paardekooper et al. (2008) have studied the dynamics and encounter velocities for 
a swarm of planetesimals in close binaries. They found inclusion of the full gas dynamics 
increased the encounter velocity and made planetesimal accretion even more difficult than in 
a static, axisymmetric gas disk. Note that their gas disk is two-dimensional, and the binary 
inclination is also not included. Nevertheless, they have shown some important results that 
we might infer the behavior in the real 3D full gas model, where both the binary 
inclination and gas evolution are included. As shown in the Figures 9-11 of Paardekooper 
et al. (2008), obtial alignments, both of eccentricity and longitude of periastron, are still 
maintained for bodies of the same size, leading very low AV between equal-sized bodies as 
in a static, circular gas disk. This means gas evolution only increase AV between different- 
sized bodies. In a real 3D full gas model, orbital alignments should be probably maintained, 
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but it is in the 3D space, and thus different-sized bodies are separated from each other, 
leading to most colhsions occurring between similar-sized bodies with low AV . Therefore, 
even in the full gas model, inclusion of a small binary inclination might be also a promising 
mechanism favoring planetesimal accretion. The only uncertainty is to what extend can it 
favor accretion, or whether it is enough to render accretion-friendly impacts dominant? This 
problem need to be addressed by a future work which can model a 3D evolving gas disk in 
an inclined binary system. 

4.2. Planetesimal Growth: Type II Runaway Mode 

The major result of this paper is: inclusion of a small binary inclination favors planetesi- 
mal accretion significantly, although it might not be enough(if only the green and blue region 
are counted as accretion) for planetesimal growth to proceed. Even if planetesimal growth 
can proceed, it cannot be classical single-star-like accretion, or normal accretion. Actually, 
normal accretion collisions never exceed 20% in all the cases we investigated. Growth of plan- 
etesimals in close binary systems is mainly through "high-A^ accretion" or type II runaway 
mode (Kortenkamp et al. 2001). Type II growth mode is characterised by an initial orderly 
growth stage followed by a runaway stage. At the beginning, is increased and held 
on values exceeding the biggest objects' escape velocities, but still leading to slow orderly 
growth. After this orderly phase, growth then switches to runaway when the biggest bod- 
ies reach a size large enough for significantly enhancing their gravitational focusing factors: 
Fg = {Vesc/AVy. This turn-off size directly depends on AV^. For example, if the runaway 
growth is turned on at a gravitational focusing factor of 1 — 100, the biggest objects' escape 
velocities should be 1- 10 times as large as AV. For the coplanar case considered in this 
paper, the average AV^ is about 200 m.s~^ at 1 AU, thus requiring a turn-off size of 200 - 
2000 km. While for the inclined case with is = l(see Fig. 3), it requires a turn-off size of 
only 50 -500 km, since the average A^ is reduced to 50 m.s~^. Therefore, inclusion of the 
binary inclination, the runaway growth is shifted to an earlier date. 

5. SUMMARY 

In this paper, assuming a small inclination of the binary orbit, we numerical investigate 
the conditions for planetesimal accretion around a Cen A. Different-sized (1-10 km randomly) 
planetesimals are embedded in the gas disk and moving around the primary from 1 to 2 AU, 
an extension of its habitable zone(~ 1-1.3 AU). The gas drag is modeled with a simplified 3-D 
axisymmetric non-evolving gas disk. Fisrt, we study the effects of the binary inclination on 
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planetesimal accretion by comparing two nominal cases with and without inchnation(zB = 
and = 1 deg). Then we explore more cases with different binary inclinations and gas 
densities. The main results are summarized as follows: 

1) In the nominal coplanar ca.se{iB = 0), planetesimal disk is almost damped to a 2D 
plane, and thus impact rate is highly increased. While in the nominal inclined case {zb = 1 
deg), bodies of different sizes are forced to orbits of different inclinations (3D-differential- 
orbital-phasing) and well separated from each other, leading to much lower total impact rate 
compared with the coplanar case. 

2) In the nominal coplanar differential orbital phasing is almost in a 2D plane, 
impacts between different- sized bodies are favored, leading to high AV and planetesimal 
erosion in all the region studied in this paper ( 1- 2 AU to the a Cen A). While in the 
inclined case(zB = 1 deg), because of the 3D-differential-orbital-phasing, impacts mainly 
occur between similar-sized bodies, thus leading to much lower AV^ and a possible accretion 
zone, which can extend to as far as about 2 AU to the a Cen A when assuming an optimistic 
Q* value. 

3) In the small inclined case(zB = 0.1 deg), the results are marginal. Different-sized 
bodies are not fully separated from each other, and thus a large part of impacts still occur 
between them, leading to conditions, although much better than that in the coplanar case, 
still accretion-hostile beyond 1 AU. While for the larger inclined ca.se{iB = 5 deg), results 
are comparable with those in the nominal case (z^ = 1 deg). 

4) For the nominal binary inclination(i5 = 1 deg), as gas density increases from 0.1 to 
10 MMSN, condition becomes less (more) accretion-hostile in most outer (inner) region. 

5) Inclusion of a small binary inclination(say 1-5 deg) significantly favors accretion 
by reducing A^. However, the balance between accretion and erosion still has a large 
uncertainty, and the reduced AV is not low enough for normal accretion. Planetesimal 
growth in close binary systems, such as a Cen A, is probably through "high-AV^" collisions 
or type II runaway mode. 

We thank Thebault, P. for discussions and valuable suggestions. This work was sup- 
ported by the National Natural Science Foundation of China (Nos. 10833001 and 10778603) 
and the National Basic Research Program of China(No.2007CB4800). 
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A. THRESHOLD VELOCITY 

The threshold velocity of the erosion/accretion limit mainly depends on the critical 
specific energy Q* . However, different Q* available in literature often significantly diverge 
from one another. As in Figure 8 of Benz & Asphaug (1999), for example, the difference in 
Q* between Benz & Asphaug (1999) and Holsapple (1994) can be as much as ~2 orders of 
magnitude in our studied size range(l-10 km). Such an broad Q* is less usable to interpret 
our results. Therefore, first, in order to have a clear erosion/limit, we adopt our nominal 
Qnom = Q*mi i'^ Housen et al.(1991), second, given the consideration of the uncertainty 
of Q* , we assume an optimistic Q^ig^ = ^Q*h9i- The Q* in the range of Q*hqi — ^Q*mi is a 
reasonable intermediate value between those in the Holsapple (1994) and Benz & Asphaug 
(1999). 

Following the detail procedure as in the appendix B of Thebault & Augereau (2007), we 
then deduce an "optimistic threshold velocity" Vugh, by adopting a critical specific energy of 
Qhigh and a hard excavation coefficient a = 10~^, and a "pessimistic threshold velocit" Viow, 
by adopting a critical specific energy of Qnom ^ weak excavation coefficient a = 4 x 10~^. 
Figure 6 shows Vhigh and Viow vs the radii of target for a projectile of radii of 1 km. 
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Fig. 1. — Snapshots of four kinds of size planetesimals(l, 2.5, 5 and 10 km) in the X-Z 
coordinate plane. Left: nominal coplanar case. Right: nominal inclined case. 
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Fig. 2. — Time-evolution of impact distribution on the R1-R2 plane. Labels of the color bars 
are the number of impact. 
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Fig. 3. — Comparing the accretion conditions for two nominal cases: is = and zb = I deg. 
The two top panels: distribution of AV, as a function of distance to the primary star. Each 
black point denotes an collision event, and the red scatter-lines are the average values. The 
two bottom panels: relative importance of different types of collision outcomes, as a function 
of distance to the primary star. The vertical lines denote the habitable zone ~ 1-1.3 AU. 
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Fig. 4. — Accretion conditions for two cases with the same gas density of 1 MMSN but 
different binary inchnations: is = 0.1 deg and is = 5 deg respectively. Left panels: impact 
distribution on the R1-R2 plane. Right panels: same to figure 3 but for the cases of = 0.1 
deg and is = 5 deg. 
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Fig. 5. — Same to figure 4 but for two cases witli different gas densities: 0.1 and 10 MMSN 
respectively. Both cases have the same binary inclination of 1 deg. 
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Fig. 6. — Threshold velocity vs radii of the target with the projectile's radii Rp = 1 km. 
Dot: a optimistic estimate of threshold velocity. Solid: a pessimistic estimate of threshold 
velocity. 
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Table 1: Initial Parameters for the Simulations. 





Paremeters 


Values 


Units 


Descriptions 




as 


23.4 


AU 


semi-major axis of companion 




eB 


0.52 




eccentricity of binary 


FP* 


is 


0, 0.1, 1, 5.0 


degree 


inclination of binary 




Ma 


1.1 


Me 


mass of primary 




Mb 


0.93 


Me 


mass of companion 


UD* 


Rp 


1 - 10 


km 


planetesimal radii 


UD* 


Qp 


1.0 - 2.0 


AU 


semi-major axes of planetesimals 


UD* 


€p 


- 10"^ 




eccentricities of planetesimals 


UD* 


Ip 


- 5 X 10-5 


radian 


inclinations of planetesimals 




Pp 


3.0 


g ■ cm~^ 


planetesimal density 


FP* 


Pgo 


1.4 X (10-1°, 10-^10-S) 


g ■ cm~^ 


initial gas density at lAU 




Pa 


^-11/4 




initial gas density radial profile 




Np 


10^ 




number of test particles 




Rinflated 


5 X lQ~\r,oi/AUY''' 


AU 


Inflated radius for collision search 




^ final 




yr 


final timescale for each run 



Note. — Mq stands for the solar mass. The fields marked by FP* are explored as free parameters in the 
simulations, while the ones marked by UD* have an uniform distribution. 



